Origin Of Current-Induced Forces In An Atomic Gold Wire: A First Principles Study 



m 
O 
O 

(N 
G 

1—5 

o 

(N 



i 

-a 
c 

o 
o 



> 
o 

o 

O 

-i— > 



i 

c 

o 
o 



X 



Mads Brandbyge, 1 'EI Kurt Stokbro, 1 Jeremy Taylor, 1 Jose-Luis Mozos, 2 and Pablo Ordejon 2 

1 Mikroelektronik Centret (MIC), Technical University of Denmark, Bldg. 345E, DK-2800 Lyngby, Denmark 
i lnstitut de Ciencia de Materials de Barcelona, CSIC Campus de la U.A.B., 08193 Bellaterra, Spain 

(Dated: February 1, 2008) 

We address the microscopic origin of the current-induced forces by analyzing results of first prin- 
ciples density functional calculations of atomic gold wires connected to two gold electrodes with 
different electrochemical potentials. We find that current induced forces are closely related to the 
chemical bonding, and arise from the rearrangement of bond charge due to the current flow. We ex- 
plain the current induced bond weakening/strengthening by introducing bond charges decomposed 
into electrode components. 



PACS numbers: 73.40.Jn, 66.30.Qa, 87.15.By 



The failure of interconnects in integrated circuits 
is often due to directionally biased diffusion of 
atoms caused by the presence of an electric current 
(electromigration) Q . In recent years the prospect 
of electronic devices operating essentially on the 
atomic/molecular scale has gained significant interest 0, 
and electromigration and current-induced conforma- 
tional changes are important issues when downscaling 
electronic components to these sizes A single atom 
wide gold contact is among the simplest atomic scale 
conductors and therefore ideal for fundamental studies 
of the current-induced forces. It is remarkable that these 
systems can sustain voltages up to several volts before 
breakingHH|||. 

Our aim here is to uncover the microscopic origin of 
the current-induced bond weakening/strengthening in an 
atomic gold wire. We analyze state-of-the-art first princi- 
ples calculations of the current-induced forces while other- 
investigations have focused on the theoretical foundation 
of these forces 0- Traditionally the microscopic origin of 
the current induced force on an atom is related to electro- 
static forces and the momentum transfer by the electron 
flow, the so called electron wind force 0. We find that 
the current-induced forces are rather linked to the re- 
distribution of bond charge or overlap population in the 
system. By extending the concept of overlap population 
to a noncquilibrium system with two electrochemical po- 
tentials, we can rationalize the first principles results in 
terms of a simple two orbital interaction model. 

We have performed our calculations within density 
functional theory (DFT]_J| @, [uj. We employ the 
TranSIESTA programfnj which allows a full atom- 
istic description of the scattering region and the elec- 
trodes. DFT has previously been used to describe the 
current-induced forces acting on single atoms 0, 
molecules Q, and an atom adsorbed on an nanotube 13j . 
Recently, the current-induced forces in atomic wires, and 
the resulting "embrittlement" , was addressed by more 
approximate theory [l^. 

We consider a simple symmetric geometry consisting 
of three atoms connecting electrodes in both the (100) 



and (111) directions, see Fig. ^ The electrode-electrode 
distance is in both cases chosen to be 9.6A and the wire 
atoms are initially relaxed at zero bias. We choose this 
particular geometry since here the middle atom (2) can 
roughly be considered to have only two bonds namely to 
the left atom (1) and the right atom (3). 

When the voltage is applied we observe a substantial 
redistribution of the electronic charge. At 0V the wires 
have a net negative charge of about 0.2 excess electrons. 
About half of this resides on atom 2. At 2V the excess 
electrons on atom 2 are almost halved. This charge re- 
distribution results in an asymmetric voltage drop where 
the main drop is between the negative electrode and atom 
2 (see Fig. HJl. 




FIG. 1: (a) Direction of forces and voltage drop for the gold 
wire connecting (100) electrodes at IV bias, (b) Isodensity 
surfaces for the change in density from 0V to IV. Dark is 
deficit and white is extra electron density. The solid (dotted) 
surface correspond to ±5 • 10~ 4 e/A 3 (±2 ■ 10 



"e/A 3 ) 



In Fig. [2K we show the calculated bond forces for the 
(100) / (111) wires with atomic positions fixed at the re- 
laxed values^} obtained for zero bias. We find that the 
forces on the symmetric atoms 1 and 3 roughly follow 
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FIG. 2: (a) The bondforce and overlap(bond) population for 
the (100) and (111) connected wires, (b) The relaxed wire 
bondlengths for finite bias. 



each other, also beyond the linear bias regime 7] . We 
note that the order of magnitude and direction of the 
bond forces at 2V are similar for (100) /(111) while the 
detailed behavior differs: the bond forces for the (100) 
wire stay linear up to IV and level off around 2V while 
we observe an onset of the bond forces for (111) at IV. 
The magnitude of the bond forces of the order of 1-1. 5nN 
for 1-2V is sizable compared to the force required to break 
a one- atom thick wire ~ 1.5nN fig ]. 

Along with the bond force we also show the overlap 
population (OP) in Fig. The OP is a measure of 
the electronic charge residing in the bond. We note that 
except for the (lll)-wire at low bias, the force and OP 
curves are almost mirrored. The result that an increase 
(decrease) in bond charge makes the bond force compres- 
sive (repulsive) is what we expect intuitively. In Fig. [2h> 
we show the bond lengths after subsequent relaxation of 
the wire atoms for finite bias. The main relaxations take 
place in the wire bonds with smaller displacements in the 
wire-electrode bonds. 

We aim at an qualitative understanding of the results 
for the dependence of the force on the voltage bias. To 
this end we focus on the OP. We can express [17] the force 
acting on atom i due to the valence electrons using the 



force operator, Fj, and the density matrix D, 

Fi = Tr[Fj D] where, F, : = ^ . (1) 

dRi 

Consider the bonding between two atoms represented by 
orbitals \4>i) and \4>2) and separated by bond length b. 
From Eq. \T\we get the bond force, 



F 



bond 



= (2) 



02 <Pl 



where we assume that only the hopping element, 
(4>2\H \4>i), changes with b. The bond force is propor- 
tional to the OP for the 12 bond, 



Oi2 — 2 S12 D12 



(3) 



The OP is typically taken as a simple measure of the 
strength of a chemical bond as suggested by Mulliken for 
molecules [lg. For extended systems the density matrix 
can be expressed in terms of the density of states (DOS) 
matrix, p. In equilibrium we have, 



O 



12 



2S 



12 



de p 12 (e)n F (s - E F ) 



(4) 



where n F and Ep is the Fermi function and energy. The 
contribution from states at different energy(e) in the ex- 
tended system to the bonding between the two orbitals 
is described qualitatively by the OP weighted DOS (OP- 
WDOS) or COOP curve, 2Si 2 p 12 (e), as discussed by 
Hoffmann[3. 

Now we consider the nonequilibrium situation where 
an electrical current is running though a contact con- 
necting two reservoirs, left and right(L,i?), with different 
chemical potentials(/iL, fin). The contact contains the 
12 bond. In this case the bond force will change since 
the density and Hamiltonian matrices deviate from equi- 
librium values. The appropriate density matrix for the 
nonequilibrium situation is constructed from scattering 
states, and the total spectral density matrix(p) is split 
into partial contributions cor resp onding to their left(p L ) 
or right(p ii ) origin (see e.g. 



rrespc 

13), 



p{e) = p L {e) + p R {e). 



(5) 



The density matrix for nonequilibrium is then found by 
filling the left and right originating states according to 
the respective chemical potentials [7J, [lfj|, 



D 



de p L (e)n F {e - fx L ) + p R (e)n F (e - fi R ) . (6) 



We split the overlap population(Eq. QJ further into its 
partial left and right components, 



Ou = O l 12 + Og , 



(7) 
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FIG. 3: The change in left and right partial OP, {Oi 2 ,Of 2 } 
and {023,023} for a change in bias from OV to IV for the 
(100) wire. These are broken into main orbital contributions 
with respect to the middle atom (2). We show the sum of the 
contributions for the degenerate d yz and d zx . 



where 



u 12 



2S 15 



de Pi 2 {e) n F (e- p L ) . 



(8) 



and likewise for R. When more orbitals are involved on 
atom 1 and 2 the overlap populations can be broken down 
in orbital components. In this case each orbital overlap 
will contribute with different weight to the bond force ac- 
cording to their different force matrix elements H'. This 



complication, together with the neglect of change in H' 
with applied voltage, makes the bond force and OP only 
roughly proportional, as we observe for the gold wires in 
Fig. |21 However, as a first approximation, we can explain 
the observed change in bond forces from the change in 
OP, which in turn means the change in left (AO L ) and 
right (AO R ) contributions. 

In Fig. |3| we consider the change from 0V to IV in 
left and right OP of the 12 and 23 bonds for the (100) 
case (we see a quite similar pattern for (111)). For both 
bonds AO L is negative (decreasing bonding) while AO R 
is positive (increasing bonding). We have resolved these 
further into the orbitals on atom 2 and display the main 
contributions in Fig. [3] 

In order to understand the physics involved we present 
in Fig. 2^, a simple picture. In the presence of current the 
electrons from the left electrode populate states above 
the original Fermi energy, Ep . These populate bonding 
or anti-bonding states depending on the position of Ep. 
Likewise there is an depopulation in electrons from the 
right. To make this picture more explicit we consider 
in Fig. 0t> a single bond between two atoms each with 
a single atomic level, Ep and er. The levels couple via 
a matrix element t < 0, and an overlap S > 0, and 
form bonding/anti-bonding orbitals, £3 and ea- The left 
and right atoms are coupled to left and right electrodes 
modelled by a level broadening of the left level and 
for the right. For this bond the resulting partial DOS 
due to the states originating from the left is, 
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(9) 



We note: (i) p R is found by exchanging L <-» R. (ii) p L 
is proportional to the coupling Tp. (iii) The cross-over 
between bonding(positive) and anti-bonding(negative) is 
given by the position of er. From (i), (iii) we see that 
pp contains more anti-bonding character than pr when 

EL > Er. 

We have plotted L/R COOP curves in Fig. |Hc-e) f° r 
different fillings and Tp = Tr. In the top row we assume 
that there is no voltage drop between the atoms, Ep = er. 
Here the change in OP comes entirely from the change 
in filling due to the shift of pp and pp, away from the 
equilibrium E F (Eq.|5J). For almost filled states (Fig.^Ji, 
top) the lowering of p r deplete states with anti-bonding 
character and O r increases. The increase in pp leads to 
a decrease in O l which does not balance the increase in 
O r . The net effect is a decrease in anti-bonding charac- 
ter and strengthening of the bond. When the states are 



almost empty (Fig. top) we also observe a strength- 
ening since now O l increases more than O r is decreas- 
ing. For a bond initially with maximal strength (Fig. ^i, 
top) and, in general, for Ef located in a region where the 
COOP decreases, the bond weakens since both O r and 
O l decrease. 

In the bottom row of Fig. 01c-e) we assume that the 
full voltage drop takes place between the two atoms, 
Ep — er, = pp — pp > 0. The bonding peak becomes 
dominating in the right COOP and likewise for the anti- 
bonding peak in the left COOP. The net result is less 
bond strengthening in (c) and (e), and more bond weak- 
ening in (d). The largest effect of a voltage drop is seen 
for half filling (Fig. ^1, bottom) where a large negative 
AOp occurs because the bonding peak decreases and at 
the same time more anti-bonding states are being occu- 
pied. 
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FIG. 4: (a) Generic picture of the COOP(x-axis) curve and 
chemical potentials of the electrodes, /j,l,Hr, deviating from 
the equilibrium, fiL = l-iR = Ef- (b) Simple model with 
a single bond between two atomic orbitals(eL, £r) forming 
bonding(es) and anti-bonding states(eA)- Atom L couples 
to the left electrode by Yl (life-time Ti/Yl) and likewise for 
atom R. (c-e) COOP curves for different initial filling factors. 
Left(Right) quantities are denoted by dotted (dashed) lines; 
zero voltage quantities are solid. The top(bottom) panels 
represent zero(full) voltage drop between the atoms, el = 
£r{£l — £r — /J-l — Hr)- 

In summary, for equal coupling, Yl = Tr, our sim- 
ple model yield: (i) bond strengthening for almost 
empty/filled states, (ii) bond weakening for about half 
filled states, and (iii) a bond with a voltage drop is weaker 
than it would be without a drop. When the bond is not 
coupled equally well to left and right, Yl ^ Tr, the left 
and right overlap contributions have to be weighted ac- 
cordingly. 

In the case of the gold wires we can to a good approxi- 
mation assume that the d yz , d zx , and p z orbital on atom 
2 only couple with the corresponding orbital on atom 1 
and 3. The d yz ,d zx states are almost filled and match 
the case in Fig. \Q: AO R > while AO L < 0. On the 
other hand the p z states are almost empty and match 
the situation in Fig.H;: AO L > while AO R < 0. This 
is exactly the behavior we observe in Fig. [3] For the s 
orbitals, which are slightly more than half-filled, we can 
expect a decrease in O l , which will be especially pro- 
nounced if a voltage drop takes place in the bond (cf. 
FigEJi, bottom). This is also what we observe for the 12 
bond in Fig. [3J 

Based on our simple model we sec that the almost com- 
pletely filled states(d) and slightly more than half-filled 



states(s) together yield a decrease in L and increase in 
O r . Since bond 12 couple more to the left electrode the 
total change in O l will be bigger than the change in O r . 
This explains the bond weakening of bond 12. The oppo- 
site is the case for bond 23. Also the voltage drop taking 
place in bond 12 will weaken this bond with respect to 
bond 23. 

In conclusion we have analyzed first principles calcula- 
tions of the current-induced forces in an atomic gold wire. 
We have shown that the current-induced forces are due 
to a shift in the population of bonding and antibonding 
levels, and can be explained by examining the change in 
electrode decomposed bond charges (overlap population) 
with applied voltage. Although we have concentrated 
on a simple system the approach should apply to atomic 
scale conductors in general. 
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